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ABSTRACT 


A complementary code pair has the property that its autocorrelation 
function contains no range sidelobes. This property makes complementary 
codes attractive for use in signal waveform design for detection and 
tracking systems. An extensive list of the properties of these codes 
has been formulated by M J. Golay and additions were made by S. Jauregul. 
However, the ambiguity functions of complementary codes have not been 
previously formulated. 

This thesis considers the anak ty functions of complementary 
codes for two cases: 

1. Narrow-band analysis 

2. Wideband analysis. 

In the narrow-band analysis it 1s assumed that the doppler effect 
does not significantly alter the envelope of the transmitted signal. 
Siebert's definition of the ambiguity function is utilized in the formu- 
lation of the narrow-band complementary code ambiguity function. 

In the wideband analysis it is shown that the doppler effect sig- 
nificantly alters the envelope of the transmitted signal. This alteration 
of the envelope of the transmitted signal becomes significant when the 
product of the two factors, the relative velocity between the transmitter 
and the target, and the time~bandwidth product, approaches one-half the 
velocity of propagation in the medium. Siebert's ambiguity function is 


extended to include the wideband case, 
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CHAPTER I 


INTRODUCTION 


The detection problem in modern search radar systems is generally 
concerned with the optimization of four major factors. These four 
factors are: 

1) Maximization of the effective range of the system. 

2) ayia of the number of false alarms due to 

impulse noise. 

3) Maximization of range and velocity resolution. 

4) Accurate measurement of the target's position and velocity 

without ambiguity. 

The solution of this optimization proglem is predicated primarily 
on a suitable choice of a transmitted signal waveform, To emphasize 
this point let us analyze the four factors that are involved in the 
optimization problem, 

The first factor can be analyzed from the standpoint of the radar 
range equation. The prominent result of this analysis is that the 
system designer may increase the range of the system by increasing the 
duration of the transmitted pulse, 

The second factor can be disposed of by setting a high enough 
threshold such that very few noise impulses exceed the threshold, but 
not so high as to degrade the probability of detection below a tolerable 
level. 

The last two factors depend almost entirely on a judicious choice 
of a signal waveform. Hence it is topical to consider a measure or cost 
function by which various waveforms may be compared, Siebert's! modifi- 


Z 


Cation of Woodward's ambiguity function provides such a comparison. 
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It is for this reason that the ambiguity function forms the basis of 
the modern radar detection philosophy. 

This thesis is primarily concerned with a class of codes, called 
complementary codes, which minimizes range ambiguities. ‘lost codes, 
which have been described and analyzed in the literature on radar wave- 
form design, have the drawback that range and doppler ambiguities exist. 
Examples of these codes are: 

ly ebinear ©) (chirp) 

2) Truncated pseudo-noise codes 

3) Barker codes. 

Of the codes mentioned above, the Barker codes' autocorrelation functions 
exhibit the smallest range sidelobes, and have been named "perfect" words. 

Complementary codes have the feature that their autocorrelation 
functions have no range sidelobes. Although much work has been done in 
the area of obtaining properties of these codes, the analysis of the 
ambiguity functions for these codes has not been done. It is the pur- 
pose of this thesis to analyze complementary codes, in both range and 
doppler, for two cases: 

1) Narrow~band analysis 

2) Wideband analysis. 

By narrow-band we mean that the doppler effect does not significantly 
alter the envelope of the transmitted signal. The narrow-band analysis 


will be formulated by utilizing Siebert's ambiguity function.* However, 





*Siebert's ambiguity function is a modification of Woodward's 
ambiguity function, The essential difference is that Siebert's ambiguity 
function is real valued, whereas Woodward's ambiguity function is complex 
valued, 


Ag 





one of the postulates in the formulation of Siebert's ambiguity function 
is that doppler does not significantly affect the envelope of the trans- 
mitted signal, Hence Siebert's ambiguity function cannot be used in the 
wideband analysis. We will extend Siebert's ambiguity function to include 
the wideband case, In Appendix B we will define the region of validity 
of the narrow-band analysis. The essence of the derivation of the region 
of validity of the narrow-band analysis is that the time-bandwidth-velocity 
product must be much less than the velocity of propagation divided by 
two, in order that we may assume that the envelope of the transmitted 
signal is not altered by the doppler effect. Then by wideband, we mean 
that doppler has a significant effect on the envelope of the transmitted 
signal. Throughout this thesis, we will assume that target acceleration 
is negligible over the interval of measurement. 

Let us now consider complementary codes. Complementary codes or 
series were originally conceived by M. J. Golay,3»4 in his work on 
infrared multi-slit spectrometry. An extensive list of properties for 
complementary codes was formulated by M. J. Golay and additions, based 
upon applications of group and vector theory, were formulated by 
S. Jauregui in his doctoral dissertation.° 

A set of complementary series may be defined as "a pair of equally 
long, finite sequences of two kinds of elements which has the property 
that the number of pairs of like elements with any one given separation 
in one series is equal to the number of pairs of unlike elements with 
the same given separation in the other series."© These sequences are 
binary and are aperiodic. 

It has been shown that the autocorrelation function of these codes 
gives rise to a central peak. The height of the peak is dependent on 


the length of the code, for zero delay, and decreases linearly to zero 


i 


for delay equal to the duration of one code symbol, then remains at zero. 
It is this property of the complementary codes that makes this coding 
method attractive for use in the detection problem. Because the defi- 
nition of the complementary codes appears rather formidable, we will 
utilize the property given above whenever we refer to complementary 
codes. Appendix F is a partial list of the properties of complementary 
codes. This list is provided so that complementary codes of almost any 
length may be generated from basic length codes called kernels. A kernel 
is defined as a basic length code which cannot be decorposed into shorter 
length codes by an inversion of the standard generating methods. To the 
author's best knowledge there are only four kernels that are known to 
exist. These kernels are listed in Appendix F. 

At this point we might indicate what the ideal ambiguity function 
would be. Since we would like to resolve two targets in range and doppler 
to an infinite degree, the ideal ambiguity function would be a delta 
function at the origin of the range, doppler axis. Thusly, the system 
would have infinite resolution. This can never be achieved; what would 
be required is a code of infinite length with the duration of each code 
symbol infinitesimally small. However, we should not discard the idea 
of an ideal ambiguity function, since it gives us a reference by which 
to compare other ambiguity functions, 

A brief synopsis of each chapter will be now given. 

In Chapter II a system model for single codes is postulated and a 
correlation between a physical waveform in the syster and the ambiguity 
function is given. By single codes we are referring to binary codes 
such as Barker codes, pseudo-noise codes, etc., i.e., codes which con- 
sist of a single sequence; as opposed to complementary codes which 


consist of a pair of equally long binary sequences. In this chapter 


14 











the analysis is entirely narrow-band. The single code amLiguity function is 
derived to ease computational matters in the derivation of the ambiguity func- 
tion for complementary codes. 

Chapter III contains the derivation and the justification of the form of 
the narrow-band complementary code ambiguity function. 

In Chapter IV a system model is postulated for the complementary codes 
and the derivation of the complementary code ambiguity function is accomplished. 
Computer programs were run for the two kernels of length 10 and the kernel of 
length 26, and ambiguity diagrams were constructed, One property of the entire 
class of binary codes in which phase-reversal modulation is utilized, is un- 
covered in this chapter. It is proved that in the doppler domain, the zero 
crossings of the ambipuity function move closer to the origin of the ranre- 
doppler axis as the length of the code is increased; assuming that the duration, 
T, of a code symbol remains constant. The implication of this phenomenon is 
that better doppler discrimination is possible for longer code lengths, all 
other things being equal, 

Chapter V contains the investigation of the doppler effect on the trans- 
mitted signal and the proof that the envelope of the transmitted sifnal is 
altered due to doppler. This chapter in conjunction with Appendices A and D 
shows that the doppler effect is threefold: 

1) The carrier frequency is shifted by the doppler frequency. 

2) The envelope of the transmitted signal is compressed or expanded 

according to whether the target is closing or opening with respect 
to the receiver. 

3) The received and transmitted energies are different; in fact, the 

received and transmitted energies differ by the same factor as the 


received and transmitted carrier frequencies. The above statement 


he 


assumes a lossless medium and ignores spherical spreading of 
the transmitted and received waves, 

Chapter VI contains the derivation of the wide-band ambiguity 
function for single codes as well as complementary codes, valid for 
positive differential doppler and positive range deviation. Unfortunately 
due to time limitations, computer programs to generate numerical values 
for the wideband complementary code ambiguity function were not acconm- 
plished. 

To clarify the discussion of the model of the system, which is 
postulated in Chapter II, a brief formulation of the autocorrelation 
function of complementary codes will be given. As mentioned previously, 
complementary codes consist of two binary sequences, which will be 
referred to as Code A and Code B or simply A and B. To obtain the auto- 
correlation function’ for the complementary code, RochiT)ag we compute 
the autocorrelation function for code A, Ry(T), and for code B, Rp(T). 


Then Ge 1s obtained by summing Ry Oe) and Re (ar) 5 ig Os, 
RC) eRe) a): 


The property of interest in the autocorrelation function of the comple- 
mentary codes is that no range sidelobes exist. All binary codes have 
autocorrelation functions with range sidelobes. However, R(T ) has 

range sidelobes which are the negative of the sidelobes of Rp(TZ ). 

Hence when R,('T) and Rp(T ) are summed, the range sidelobes cancel. Fig- 


ure 1-1 shows this unique property of the complementary codes. 
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Fig. i=-1a - Autocorrelation function for code A, length 10; 
A = set £-1,41,4+1,-1,4+1,-1,+1,+1,+1,-1} 


Fig. 1+1b - Autocorrelation function for code B, length 10; 
Bo = Sub 4-1 gabilgatatvyeetsl , tke, ale, -F 1) oeateed ata 
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CHAPTER II 
MATCHED FILTER APPROACH FOR THE 
SINGLE CODE SYSTEM (NARROW-BAND ANALYSIS) 

It might be well to look at the single code system before attempting 
to analyze the complementary coded system. In the matched filter 
approach the underlying rere is to correlate the ambiguity function 
with a physical waveform in the system. The ambiguity function is an 
excellent measure of range and velocity resolvability, 1.e., how closely 
one can distinguish two targets in range and velocity. A common mis- 
conception is to state that the output of the matched filter is identical 
to the ambiguity function for all cases. This statement holds some 


truth; however, it is not generally correct in carrier modulated systems. 


Pe ~*~ 


Let us depart from this discussion to model our system before attempting 
to clarify the above statement. The model that is chosen is a slight 


modification of the model used by Davenport. ° 


[Transmitter | s(t). Transmitted | s'(t-td) + n(t) _| Receiver 
, Medium | = 
i re 


t 


| information’ idisturbance| 
source | (additive) 
(target) | 


System model 


Figure 2-1 


The difference between Davenport's model and the model that is chosen 
is that in the radar, as well as the sonar case, the information source 
interacts in the medium, whereas in Davenport's system model, the infor- 


mation source interacts in the transmitter, 
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The transmitted signal is chosen as a phase-reversal modulated 


signal, where 


S(tjes) Bx coaGum et — nt ) 2-1 
where ine = j0 depending on the value of the ith code 
| 
; symbol 
\ = normalization constant. 


For convenience the signal energy is normalized to unity by specifying 
that 


\s*(t) at = 1 : 


It is a well known fact that a phase-reversal modulated signal is 


identical to a DSB/SC signal; hence the transmitted signal can be written 


as 
si(t) = Ct) cosmo t, 2-2 
where L(t) = MeN, 17 (tein) - ie - G+ at, 
ra 4.0 
and Xj; = - 1, depending on the value of the ith code 


symbol, 


The normalization constant can be easily evaluated by squaring equation 


2-2 and performing the integration. Thus 





. 2 NI 
ee 2 : . 
\s2(t) dt ar \ >t at eeaere.t.,| dt 2-3a 
ae : 
ean \e , Agin2 w oNT yecie 
2 4 G9 ° 


The second term in equation 2-3b is very small compared to the first 


term, for ., large. This term can be interpreted as the difference 


re 


between the initial and final value of the carrier phase and will be 


small, Hence the normalization constant is approximately given by 


———_ 


2 
Oo “|e 2=3¢c 


4 


The received signal contains a noise term, n(t), which will be 


9 


ignored since the receiver will utilize a matched filter” to raximize 
the signal-to-noise ratio. Thus the received signal, 5 (eae. is “tie 
transmitted signal modified by the information source, and contains 
range and velocity information. Let us assume that the target speeds 
are sufficiently less than the speed of propagation of the wave in the 
medium (say of the order of 1074 x c, where c is the velocity of propa- 
gation.) The principal doppler effect, in this case, is to shift the 
carrier frequency by the amount fy where fy is the doppler frequency; 
doppler will have little effect on the duration of a code symbol, T. 
(See Chapter IV for the discussion of the doppler effect on the curation 
of a code symbol.) With the above conditions postulated, the receivec 


Signal may be given bv 


S<+-t4) = 1 (t-t3) COS [Ku 4 wd; XA-ty A 


i 


N 
OD Xf ty T)- Atty vT] cosflas pity, 2-4 
1-0 4 


where the convention is chosen so that the doppler frequency is positive 


for closing targets, and negative for opening targets. 


At this point the question might be raised as to what are the 
various waveforms in the receiver. As mentioned earlier in this chapter, 
a matched filter will be utilized in the system model. To obtain the 
signals to correlate with the received signal, we translate the trans- 


mitted signal by various expected doppler frequencies as indicated in 
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Figure 2-2a Receiver block diagram 


P-2<b ith doppler rilter of the filter bark 
in Figure 2-2a 


a 


Figure 2-2a, The doppler shifted signal is then delayed. The delay is 
required to account for the propagation time to and from the target. To 
instrument this system in real time, the delay mechanism might be of the 
form of a tapped delay line. We shall denote the delayed, doppler shifted 


Signal as the replica signal. The replica signal is given by 
si(t - T,) = I(t - Ty)cos [( Wo + Way) (t - TEI 2-5 


Let us now analyze the various signals in the receiver. To simplify 
the analysis, the matched filter block diagram is duplicated from 


Figure 2=-2b. 


The output of the matcled filter is 


YyeC%e 54 Way) = ( S'(4-t3) $;(4-%) at 2a 


-Cco 


if 


(rats Cos (C0, + Wiens 4) | , 


TU-%) eos (wd, Wai 4 "G)] At 2-6b 

Then let 
Ee ed 2-78 
Dat Oo ge 2-96 


Equation 2-6b becomes 


O 
ad — 
Up Tn Myris) = STADIA, $4 - Tr) COS (td,,- Wa; Hy dts 


Cos | Welt ttq = T,)1 dt, 2-8 


Ze 


Dp asuel 42H pu 
O Uae kom 


rez oungty 


ISAT LF 





Dayo Tey 


Then let 


= —_ 2a 
CT ce ty 9a 
io = Oy -— Wd 2.9b 


If the above changes in variables are effected, then 


Yy CT Wy) = \ Tt) 1-7) COS (Wo, + W4)t, COS w,, Coy- T) Gb, 2-10 


The variables T and Wa deserve some comment. From equation 2-9a, 
we see that if T is zero, then the propagation delay to the target 
has been matched in range bin The Hence % is the range to the 
target. From equation 2-9b we see that if Wg is zero, then Wg, is 
the doppler frequency. Thus both 7 and Wg represent the departure 
of the target's range and velocity from the estimated range and velocity. 
The variable t, is just the variable of integration, and since it is 
inconsequential as to what nomenclature is ascribed to it, we shall 
denote ty, by t. Ww ,4 represents the estimated carrier frequency of the 
replica signal. Hence the output of the matched filter in the syster 
1s given by 

oo 


( LG) eee) cos ( Wd, + Wa)t Cos Udo, (t- T)dt : 2-lla 
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The second term in equation 2-llb has frequencies centered about twice 
the carrier frequency plus the doppler frequency, and will be rejected. 
Note that the output of each matched eee an ordinate in the T 
domain. 

Next we require a device to collect the output of each of the 
matched filters in an orderly fashion. Hence the use of the tern, 
collator, to describe that device is appropriate. We shall define the 
collator in Figure 2-2b as that device which collects the various out- 
puts of the matched filter sequentially as a function of range deviation, 
"". To clarify the above definition, we might think of the collator as 
a graph plotter, where ordinates (outputs of the matched filters) are 
plotted as a function of range deviation, ~. The "envelope detector" 
in Figure 2-2b is not an envelope detector in the normal sense, and 
requires some explanation. As mentioned above, the collator collects 
the outputs of each matched filter and produces a continuous curve in 
the T domain. This curve oscillates at Wo, in the TC domain, but only 
its envelope has information content. The "envelope detector" refers 
to that device that takes the envelope of the output of the collator, 
but in the range deviation or T domain. In the usual sense of the 
term, envelope detection is accomplished in the time domain and hence 
a distinction must be made in our usage of the term envelope detector, 


Now let us look at the output of the collator; 


on 


a 
20 


If the Euler identity for the cosine function is utilized, the output 


from the collator is 
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44.7 + st 
(Ty Wy) = ee \ Tis)T-ae' * dt 
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If we note that the second term of equation 2-13 1s the complex conjufate 
of the first term and that the sum of a quantity, say z, and its corplex 
conjugate is 


Z+ Z* = 2 Re(Z), 


then 


a Wo; 'T 


1 “y { peas { 
yCe, ary) = i Ke. 1% T4)TIG-O¢ dt. 2-14 
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The terms I(t) and I(t - 7% ) are the modulation functions and involve 
the limits of integration and the values of the codes, The modulation 


functions, according to equation 2-2 are 


N 
TH) = X2 %{f1ear-fe-eork, —, 2-15a 
N 
TET) = & 7, %; 4 lt-e-i7) — W((¢-v-G4T) { 2-15b 
ied ‘ 


The unit step functions in equations 2-15a and 2-15b will determine =i 
limits of integration and the noes are ae depending on the value of the 
codes. Figure 2-4 is the plot of the time waveforms of I(t) and I(t —"e 
To obtain the limits of integration, we merely observe the zero crossings 


of the functions. Hence the output of the collator is given by 
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The above equation can be generalized to a code of any length NM, and 
will be written as the sum of two series, Therefore the output of the 


collator is 


7 Jj Wert ( a i“ 
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where N length of the code 
k = smallest integer > | c/-- \ 


T = the duration of a code symbol = clock period 


The first integral in equation 2-16 can be evaluated as 
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If the following trigonometric substitutions are utilized, 





GOS X- GOS y = - 2 sin (AZ) sin (4) 


os 


miei ys 0) = ¥ - rx i 
«. sh\ — ee ce | n (3. ) St . { fs ) 


” 


equation 2-17b becomes 
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The second integral in equation 2-16 can be evaluated as 
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We then substitute equations 2-18b and 2-19b into equation 2-16. 
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If the exponential terms are combined, and the real part of equation 


2-20 taken, then 


2 Nk 
_ & 2 uy (Oc+4)T-T) . 
Y (7, wy) — Wy XX, Sin a 


f=4 
we N-(+1) 


os fanraGarape] ah OPE 


j=4 


ary WdCC-KT) poe [ualassoT + (2ud+ Boe | . 
A 2. , 


In order to plot equation 2-21, it is necessary to obtain the equation 


in the form of an envelope term multiplied by an RF term. Let 
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We now utilized the trigonometric identity 


cos(a + b) = cos a cos b = sin a sinb, 


Equation 2-23a becomes 
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Thus the output of the collator becomes 
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where %@ is the phase due to the amplitude of the sine and cosine terms. 


If equations 2-22 a, b, c, d, e, are substituted into equation 


2-23c, and if the sine terms are placed in the familiar sin & form, 
x 
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Equation 2-24 is plotted in Figure 2-5; note the RF fine structure 
under the envelope. 
Now let us analyze the ambiguity function. Essentially, definitions 
and results will be stated from Siebert's paper+ with the exception of 
a few terms which will be modified to conform to the nomenclature of 


this thesis. The ambiguity function as defined by Siebert is 
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where S, (t) is the complex envelope of the transmitted signal. 


In terms of this thesis 
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I(t) is real. Therefore 
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Equation 2-26 departs from equation 2-14 in that the carrier phase term 
is absent and the magnitude of the integral is taken, Utilizing the 
same analysis used in obtaining the limits of integration for the output 
of the collator, we see that equation 2-26 can be written as 
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If the results of equations 2-17 through 2-20 are utilized, equation 


2-27 becomes 
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The term te can be factored from equation 2-28, Note that 








If the magnitude of equation 2-28 is taken, then the equation in its 


final form is 
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In the earlier portions of this chapter, the statement was made 
that the output of the matched filter is identical] to the ambiguity 
function. We are now prepared to discuss this particular point. In 
comparing equation 2-24 (the output of the collator) with equation 2-29 
(the ambiguity function), we note that the ambiguity function is identical 
to the envelope of the output of the collator. Hence the difference 
between the ambiguity function and the output of the matched filter is 
that there exists an RF fine structure of frequency - > in the output 
of the matched filter which is absent in the ambiguity function. ‘ote 
that if the phase of the RF fine structure is some odd multiple of > 
then the output of the matched filter will be zero. Thus, we can say 
qualitatively that the ambiguity function measures the positive maximums 
of the output of the matched filter per cycle of the RF fine structure. 
To circumvent the problem of outputs of the matched filter being zero 
near the central peak, we sum and take the envelope of the output of 
the matched filters. However, care must be taken in choosing the differ- 
ential range deviations. Suppose that the differential range deviation, 


AT , 18 some multiple of the inverse of twice the frequency of the RF 


fine structure, 1.e., 
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Furthermore, suppose that in any matched filter, the phase difference 
between the transmitted signal and the received signal is an odd multiple 
of ae radians, Then the output from all the matched filters will be 
zero, This point 1S easy to see because the system will be tracking the 
zero crossings of the RF fine structure. 

The RF fine structure 1s also troublesome in the area of multi- 
plexing. This point will be discussed in Chapter IV, when we will 
formulate ideas on the type of multiplexing to be used in the complemen- 
tary coded system, 

The output of the “envelope detector" is given by the envelope of 
equation 2-24, However, as mentioned earlier, the output of the "enve- 


lope detector’ is identical to the ambiguity function. 


B/ 


CHAPTER Iii 
DEFINITION OF THE AMBIGUITY FUNCTIO}! 
FOR COMPLEMENTARY CODES 

In this short chapter, we shall extend Siebert's Ambiguity function 
to the ambiguity function for complerentary codes. Let us analyze the 
signal processing techniques required to produce the desired results in 
the complementary coded system. In Chapter II, it was shown that in 
the single code system, the output of the “envelope detector" in Figure 
2-2b was identical to the ambiguity function. What then, should be the 
configuration of the complementary coded system? Should the complementary 
coded system be, in essence, two single coded systers such that the 
outputs of the two "envelope detectors" be summed together? Or should 
the two codes be summed elsewhere? 

Let us postulate the complementary coded system to be comprised of 
two identical single code systems with the outputs of the "envelope 
dectors" of code A and B summed together. In this case the ambiguity 
function of the complementary coded system will be the sum of the ambi- 


guity functions of code A and BP, i.e., 


Cr, W,) ae Bae ws) ee 04) ; Jem) 


One of the properties of the ambiguity function is that 
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Heteemcaiati0n 35l.f0r Daa 0 becomes 


B.cqapl = ARM) 4 oll RC . 3-2 


= 
As indicated in the introduction of this thesis, the desired property 
of the complementary codes was that off the central peak, the autocorrelation 
function of code A, Ra(T), was the negative of the autocorrelation function 
of code B, Ry (T )s so that when the autocorrelation functions were summed, 
the range sidelobes canceled. Since the proposed definition (equation 
3-1) of the ambiguity function contains two positive quantities Ira Cz )\ ; 
and | race )| , the definition is inadequate. : 


As a second approach to the problem, let us define the ambiguity 


function of the complementary coded system as 


| -° wt 
B.(re004) = + \\ TaWIye-ve’ dt 
WAL 
* (reurtett-2 a” k 3-4 


where I,(t) = envelope of code A 


Ip(t) = envelope of code B, 


For zero doppler frequency, equation 3-3 becomes, 


¥ Cr, 0) = Of, \iR ec) +Ro(e)| : 3-4 


3 


It is easy to see that the range sidelobes of code A will cancel the 
range sidelobes of code B so that the range sidelobes for the comple- 
mentary coded ambiguity function will be nonexistent for zero doppler 
frequency. Equation 3-3 will be chosen as the definition of the ambip- 
uity function for the complerentary codes, 

Equation 3-3 gives us the clue as to the confiruration of the 
complementary coded system, We can infer that the two codes should be 


summed at the output of the matched filters, then envelope detected. 
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CHAPTER IV 
MATCHED FILTER APPROACH FOR THE COMPLEMENTARY 
CODED SYSTEM (NARROW=BAND ANALYSIS) 
A. BACKGROUND 
In the complementary coded system, some form of multiplexing is 
required since two codes will be utilized to obtain the desired cancel- 
lation of the range sidelobes. Of the three obvious types of multi- 
plexing (time division, frequency division, and quaternary coding) only 
time division, TDM, and frequency division multiplexing, FDM, will be 
discussed. 
Bi. ROPOSED FDM SYSTEM 


Assume for the moment that the FDM system is modelled as below. 
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Complementary coded FDM 
transmitter block diagram 


Figure 4=] 


It might be well to observe the transmitted spectrum to aid in the 


understanding of the signal processing schemes that will be required in the 


receiver. The frequencies, f, and fos are at RF. The outputs of each 


of the phase-reversal modulators can be easily obtained from equation 2-2 
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It is easy to see that the spectrum of the spectrum of the transmitted 
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signal consists of two functions centered at f, and fo. 
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One-sided spectrum of the transmitted signal 


Figure 4-2 


Thus to detect each code, the system requires bandpass filters centered 
at the carrier frequency plus the doppler frequency. Hence the syster 
block diagram is given by Figures 4-3a and 4-3b. Since the output of 
each range bank can be treated the same as in the single code syster, 
equation 2-15 applies with the substitution of appropriate subscripts. 
Let ®% be the envelope quantity in equation 2-15; hence the output from 


any range bank in Figure 4-3b for code A is given by 
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ca id) = ie cos| ( Wy ke ax) 'C dy, \ : 4-3a 
The output from the corresponding range bank for code B is 


Yn erty) = So cos\ (10, 4 at) TT ~ , \ P 4—3b 


The output from the summer in Figure 4-3a is the sum of the outputs 


from the range banks of code A and B. Therefore 
4% wa) = U bts wa) + Ue Ct, wd) 
= Ba cos (w+ 2h)t~ dy 


+ %e cosh (wW,+ M4 v- dal . 4-4 


The effect of the different frequencies of the RF fine structure for 
code A and B is to lower the output from the matched filter at the 
central peak and to increase the range sidelobes, It is true that, if 
the carrier frequencies of both codes are nearly equal, the range side- 
lobes will be very small; and in fact, might be quite tolerable. lowever, 
bear in mind that the characteristic of interest in a complementary 
coded system is for range sidelobes to be nonexistent for zero doppler. 
In this system the factor which is troublesome is the RF fine 
structure which is absent in the ambiguity function. Hence the analysis 
provides a further contradiction to the statement that the output of the 
matched filter is identical to the ambiguity function. Indeed, it is 
that factor which is different between the output of the matched filter 


and the ambiguity function which defeats this method of multiplexing. 


Cee TDNeSY STE 

Another multiplexing method is TDM. The suceeding portions of 
this thesis will dwell entirely on phase-reversal TDM systems for comple- 
mentary codes. The results from Chapter II concerning single code 
systems can be easily extended to the complementary coded systers. In 
the transmitter the essential difference is the time multiplexing circuitry; | 


hence the transmitter can be modelled as below. 
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TDM transmitter block diagram 


Figure 4-4 


As in Chapter II the transmitted energy will be normalized to unity 


by specifying that 
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To obtain the normalization constant, we note that 
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Therefore the normalization constant, (Q,, 1s approximately 
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In the transmitter, code B is delayed and summed with code A, The 
receiver is modelled in Figure 4-5. To detect the signal the gating 


circuit in the receiver separates the two codes. The output from each 


47 


requency | 
oft) 


Sault) “ | Saad) 
de lay] S(t-tu¥ range r Yale Wd) 


line bank 
Ss lt=tdu) range Yale wai) 
; bank [7 
Gating | 
ear Cite, ranfa Ypa(C,da) 
Dank 


summer Enve lope KCC, Win) 
Detecto 


X(€,942) 












9e(t-td) [range 
; bank X9(t,wy| Summer 


range 
bank 
range 
bank 


$3, + SG) e Senlt) 


Frequency / 


Envelope 
Detecto 


rnvelove 
Detecto 


X (eds) 









Ye (7,wd:) 








46(7%,wd,) 


to display Comparator 
Lae Ss 


TDM receiver block diagram 


Figure 4-5 
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correlator in Figure 4-5 is identical to the output from a correlator 
in the single code system. Hence the output from the “envelope detector" 


is piven by 


oe a W0;) = Envelope of 4 y,tey Ww) Sy Yat% wp} 2 Laas? 


Upon examination of equation 2-15, we find that the only terms which 
depend upon the value of the codes are the X.'S. The rest of the terms 
are common to both ya (7, Wq) and yp (7, Wg). Hence if the signed 
products in equation 2-15 are replaced by signed products for code A and 


by Ls@e, 


X;X; is replaced by (XA;XA 


; jet XByXB,) 


and if the envelope of equation 2-15 is taken, then the output of the 


"envelope detector’ will be given by 
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As stated in Chapter II the output from the "envelope detector" is iden- 
tical to the ambiguity function. Therefore the ambiguity function for 
the complementary codes will be given by equation 4-9. 

Let us now analyze the narrow-band ambiguity function for comple- 


mentary codes. For zero doppler equation 4-9 reduces to 
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Thus we see that from the analysis done in the introduction, the ambiguity 


function for zero doppler contains no range sidelobes. We note that 


this property is a direct result of the RF phase synchronization assumed 


in equations 4-6a, b. 
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For the case of zero range delay, the ambiguity function reduces to 
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Then if the sine terms are factored, equation 4-11 reduces to 
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If the Euler identity is utilized 
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Thus equation 4-12 can be written as 


A 


‘Elo, uy) a X, 











N N 
a a Ce Cos iT + \Z, sin sit] 4-13 
2. 8 


If the following trigonometric formulas?? are utilized 


N 
ee Csio = cog WAL cseom 
gb Zz ZA 
“ (Nts 
2s sin de 3 in SOE eosee - 
and 
N N 
pa Cosi + \ 2s Sint@ = gin N& cosec ° 
Yee g24 Zs hh 


eee Qian r sin (N41 | 
Zz = 


ae > 
thus equation 4-13 becomes 
» Na) QUT 
4 ae 
$.(0,)) = Of + Sn ET HT ssn Nail cos ue A 4-14 
sur ss e 
A 


Then noting that the absolute value of a product is the product of their 


absolute values, and 


je | =. 


Finally we see that equation 4-14 becomes 


B.Co,a4) = We NT ain SST 





Nyt | - 
2 


2 





If equation 4-15 is plotted, one interesting result is observed. 





Ambiguity function for complementary codes, 


in the doppler domain, 


Figure 4-6 


The first zero crossing of the ambiguity function occurs eee succeeding 
zero crossings occur at multiples of ae Hence the coded pulse sequence 
is acting like a single pulse of length NT in the doppler domain. For 

a better insight into this phenomena, let us review the origin of the 


sin MOd term. A glance at equation 4-12a shows us that the term origi- 


nated as a sum of exponentials, i.e., 


N a N {uit 
eV = i 2c A216 
gah ae e 


The criterion for the ambiguity function to be zero in the doppler domain 
1s 
A 


Sin Nay? — =, Nw Si = AyD, ene 4-AJa 
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Or alternatively, 


’ e 


Noy = yr : =o Zee 4-17b 


Hence we see from equation 4-16 that when the total phase shift is a 
multiple of 2W , the ambiguity function for zero range delay will be 
zero. This effect can be viewed as the addition of phasors such that 
when the function is zero, the phasors have added to form a closed loop. 
This, of course, occurs when the total phase shift is some multiple of 
21 . 

One added comment should be made concerning the zero crossings of 
the ambiguity function in the doppler domain. Since the zero crossings 
occur at ar, then, as the length of the code increases, the ambiguity 
function will be compressed closer to the origin; hence approaching the 
ideal. However, the length of the code cannot be increased indefinitely 
without complications. It will be shown in Chapter VI that longer codes 
will be more sensitive to ciate caused by the doppler effect, and 
unfortunately the length of the code cannot be increased indefinitely 
without reducing the output from the matched filter. The above compression 
effect on the ambiguity function is not a property possessed solely by 
the complementary codes. For example, single phase-reversal codes have 
ambiguity functions which experience the compression effect for long 
codes. It is the length of these codes, rather than their internal 
structure, which produces compression in the doppler domain; assuming 


that the duration of a code symbol, T, is held constant. 
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D. SUMMARY OF RESULTS 

Equation 4-9 was coded into Fortran and programs were run for the 
two kernels* of length 10 and the kernel of length 26. Plots of the 
ambiguity function were made for constant "doppler phase,’ ®, , where 
this is defined as the recurrent argument in our equations for the ambi- 


guity function, l.e., 


Gy= SE = AT , 4-18 
where T = duration of a code symbol. 
To convert doppler phase into radial velocity, we recall that for 


the narrow-band case 





Therefore, 


Cc f, C §, 
VY.=7- =~ OCU a 4—]! 
ar ei. Dir, = 
Table 1 gives the conversion from doppler phase to radial velocity for 
various values of carrier frequency and duration of a code symbol. 


Doppler phase, Oy » was defined in equation 4-18, so that the ambiguity 


function need not be recomputed for different values of o and T. 





*As stated in the introduction, a kernel is a basic length code 
which cannot be decomposed into shorter length codes by an inversion of 
the standard generating methods. Conversely, complementary codes of 
longer lengths can be generated by appropriately combining the kernels. 


Do) 


Figure 4-7 is the ambiguity diagram for the kernel, A, and Bo, 
code length 10, (The code values for the kernels are listed at the end 
of Appendix F.) The ambiguity diagram is the term used to describe the 
three dimensional figure generated from the ambiguity function. The 
bandwidth* of the ambiguity function in Figure 4-7 is approximately 
O.lw. For the system postulated in Figure 2-2, this bandwidth corresponds 
to a velocity discrimination of 3000 knots, for T = 1 psec, f, = 10 Giz: 
For f£, = 1 GHz and T = 0.1 psec., velocity discrimination is roughly 
300,000 knots. Hence complementary codes of this length cannot be used 
for situations requiring good velocity discrimination. 

Figure 4-8 is the ambiguity diagram for the kernel, Ap, and Bo, 
code length 10, and with @, in increments of 0.125T. 

Figure 4-9 is the ambiguity diagram for the kernel, A, and B,, 
code length 10. Note that the range sidelobes extend out to seven units, 
whereas the range sidelobes in Figure 4-8 extend to 10 units. However, 
the amplitudes of the sidelobes for the kernel, A, and By are slightly 
larger than for sequences Ay and Bo - 

Figure 4-10 is the ambiguity diagram for the kernel, A and B, length 
26. The range sidelobes extend to 14 units. The maximum amplitude of 
the range sidelobes is 0.32. The ambiguity diagrams in Figure 4-9 and 
4-10 are in the same increments of Oy » Note that the bandwidth of 
the ambiguity function for the kernel, A and B, length 26, is less than 
the bandwidth of the ambiguity function for the kernel, A, and B, length 
10. This effect is the compression phenomena for longer codes which is 


discussed in Chapter IV. 


*The bandwidth of the ambiguity function is here defined as doppler 
phase (converted to frequency) at the first zero crossing of the ambiguity 
function for T = 0. 
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TABLE 1 


CONVERSION FROM DOPPLER PHASE, @4 , TO 


TARGET'S RADIAL COMPONENT OF VELOCITY, Vrs 


Note: c = velocity of propagation of the wave in the medium 
(5.91 x 10° knots) 
ee carrier frequency 
T = duration of a code symhol 
fy = doppler frequency 
Ti = Wal = 2vyf5tw/c 








A Re OT OG LA ON OI A EN Ns ES mE es OS Ol ne 
* 


Mi nats) fe (GHz) 7 Cp ae) ; O q ar rad)| 


s 
TE LLL AEE PP TPCT 2 a. EE RE TR te I EPO TH EO vant 












t -4 
; 1 1 2.00x10 
3.38x107% 
| 2.00x1073 
3.38x1072 
1600 540x107 
2000 | 6.75x1073 
3000 | | | hes 00x1072 
1, 35x10, 
| 1.69x10 
¥ 2.00x1072 
0.2 - 
ee cee - ti nem 
| 1 2.00x1079 
| 100 | | 3.38x1075 
600 : 2.00x1074 
| 1000 | 3.38x1074 
1600 5.40x1074 
2000 6.75x1074 
3000 1.00x1073 
| 4000 1.35x1073 
5000 1.69x1073 
6000 2.00x1073 
60000 | 2.00x1072 
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TABLE 1 








(cont'd) 

60 10 i 2.00x1073 
100 3.38x10-3 
600 2.00x1072 

1000 | 3.38x1072 
1600 5.40x1072 
2000 | 6.75x1072 
3000 0.1 

4000 0.135 
5000 0.169 
6000 0.2 
60000 2.0 


60 «| 10 O.1 2.00x1074 
100 3.38x1074 
600 2.001072 

1000 3.38x107-3 
1600 5.40x1073 
2000 6.75x1073 
3000 1.00x1072 
4000 1. 3Seaio=2 
5000 1.69x1072 
6000 2.00x1072 
60000 032 


— we a 
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CHAPTER V 
INVESTIGATION OF THE DOPPLER EFFECT 
ON THE TRANSMITTED SIGNAL (WIDEBAND ANALYSTS) 


A, BACKGROUND 


In Chapters II, III, and IV the analysis was done on pulse-coded 
systems for the narrow-band case. What is meant by narrow-band is that 
doppler has negligible effect on the envelope of the transmitted signal, 
but shifts the carrier frequency by the doppler frequency. It is the 
purpose of this chapter to discuss the doppler effect on the envelope 
of the waveform and to derive the expression for the time waveform of 
the received signal (wideband case). As stated in the introduction, 
target acceleration is assumed to be negligible over the interval of 
measurement. 

B. EFFECT OF DOPPLER ON THE TRANSMITTED SIGNAL'S FREQUENCY COMPONENTS 

Assume that s(t) is any arbitrary transmitted signal. Then the 
Fourier transform (if it exists) is a measure of the frequency content 
of the signal. It is assumed that s(t) is sufficiently well behaved 
to have a Fourier transform. If the shorthand notation is utilized, the 


Fourier transform of s(t) is 


FLslay] = O(f) . 


It is well known that the doppler effect is to shift the frequency 
components of the signal such that the new frequency f' is given by 


a { . a-1 


LA4N/e 


en TLR SI NT LA 


*This equation is a statement of the two-way doppler effect for 
radar and sonar. For the radar case Torres!! has indicated that if the 
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where v_ = 
i 


radial component of velocity between the target and the 


transmitter 


C speed of propagation of the wave in the medium 


Note: The convention is chosen so that v, is positive for closing 


targets and negative for opening targets. 


Define a quantity, the doppler factor, D, such that 


Pane 5-2a 
Alternatively, 
p =_i* Wie 5-2b 
a= NC emee 


The doppler factor may be given in terms of the doppler frequency fa 
Note that the doppler frequency is defined as the change in the carrier 


frequency due to the doppler effect. Hence 
fq = i =i, 5=-3b 


Then 


fmt fo see ig Jae 





velocity vectors of the transmitter and the target are noncollinear, the 
angle at which the energy arrives at the receiver (assuming that the same 
antenna is used for transmitting and receiving) differs from the angle 


at which the energy was transmitted due to the angular aberration phe- 
nomenon. 1 


In the sonar case equation 5-1 applies if the medium is at rest 
with respect to the transmitter. 
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The above equation may be rearranged so that 


p=1+4d 5-4a 
o 


nv 
If the ratio = « 1, then the doppler factor may be approximated by 


p2fa¢ ae 5-4h 


Hence the doppler frequency is approximately 


. 4yr : 
f 4? ete 5-5 


Next let us examine the received signal's Fourier transform, S'(f). 
(The prime notation is used to indicate that the received signal differs 
from the transmitted signal.) Since the effect of doppler is to shift 
all the frequency components of the original signal by a multiplicative 


factor, D, the received signal's transform is given by 
; * 
S’ (Cf) -= SCD) ee 5-6 


C. TIME WAVEFORM OF THE DOPPLER SHIFTED SIGNAL 

It was established in the previous section that the received signal's 
Fourier transform was related to the transmitted signal's transform by 
a change in argument. Hence the inverse transform of the received sig- 


nal's Fourier transform will yield the time waveform of the received signal. 





*We arrived at the amplitude factor being unity from considerations 
of the special theory of relativity. In Appendix A we proved that the 
transmitted and received energy is dependent on the relative velocity 
between the transmitter and the target. This result was also obtained 
by A. W. Rihaczek. 12 To obtain the energy of the received signal, we 
integrate the squared magnitude of equation 5-6 over all frequencies, 
Thus the energy amplitude factor is introduced via the change in fre- 
quency scale, 
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Since the time waveform and its Fourier transform form a one to one 
correspondence, the time waveform obtained by the inverse transform of 
the received signal's Fourier transform will be unique. If the inverse 


Fourier transform is taken, then 


s(t) = Tr [8($)| = ( schyeh as ; se 


By a change in variable 


d= 


then equation 5-7 becomes 


a(t) = of Sse saa 


oh a 5-8 


Equation 5-8 has the same form as the transmitted signal with the excep- 
tion of the amplitude factor D and the change in the time variable to 


Dt. Therefore 
~{ 
s(t)= Ds(ot)= F | sG)\ 5-9 


Equation 5-9 is a well known transform pair and is referred to by 
R. MN. Bracewell ee the "Similarity Theorem." Since the received sig- 


nal is of the same functional form as the transmitted signal except for 
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a change in argument and in amplitude, s'(t) may be written as 


s(t) = DI(Ds) ces W,D4 5-10a 
N 
= Da 2 X {fio > eee wrt cos W,Dt , 5-10b 


Note that the unit step function A(dt-i7) may be written as filoi4-i 5) and 


Alote-ipil= Me-s.5) , 


Define 
I= T/D, al 
Then 
WD = wo(t + 2) = Deady . 5-12 
Hence equation 5-10b becomes 


NU 
cs! (+) = Do 2X4 | fl (4-iT) - {t-te} Cos (Wd. Oy) 2 5-13 


Equation 5-13 is the time waveform for the wideband doppler shifted 
signal. Note that the wideband signal reduces to the narrow-band signal 


if the doppler factor is nearly unity. 
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If equation 5-13 is examined carefully, we see that doppler produces 
three effects. 

1) The carrier frequency is shifted by the doppler frequency. 

2) The duration of a code symbol is altered by the doppler factor. 

3) The received energy differs from the transmitted energy by the 

doppler factor. 

The reason for the carrier frequency being shifted by doppler is obvious, 
and Appendix A shows the relationship between the received and transmitted 
energy. However, why is the duration of a code symbol altered by the 
doppler effect? 

Let us determine the physical significance of the alteration of the 


duration of a code symbol. 


Let 
nl = number of cycles of RF in the duration of one code symbol 
of the transmitted signal 
q' = number of cycles of RF in the duration of one code symbol 
of the received signal, 
Then 


q = duration of one code symbol for the transmitted signal 


duration of one cycle of the transmitted signal 5-l4a 


= a = f T 5-14b 


ie . duration of one code symbol for the received signal 5-l5a 


duration of one cycle of the received signal 


_ as 1S 
vere oO TD. , 


Oo 





From equation 5-11 
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and from equation 5=3b 


oO o* 
Therefore 
: J 
A=" (Eo) G) = fol « 5-16 
Hence 
j 
4| =" ' 5-17 


Thus we see that the number of cycles of RF in the transmitted and 
received signals remains constant. Hence the implication of the doppler 
effect is that when the transmitted signal's carrier undergoes a shift 
due to doppler, the duration of a code symbol changes so that the total 
number of cycles of RF for the transmitted and received signals remains 
constant. This is compatible with the special theory of relativity as 


shown in Appendix D,. 
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CHAR Ts val 
WIDEBAND AMBIGUITY FUNCTION FOR 
COMPLEMENTARY CODES 


A. BACKGROUND 


It might be well to indicate those situations in which the doppler 
effect on the envelope of the transmitted signal cannot be ignored. 
Certainly for relative velocities, between the target and the transmitter, 
near the velocity of propagation of the wave in the medium, this effect 
cannot be ignored. Another situation involves the duration of the pulse 
train. As indicated in Chapter V, the doppler shifted waveform undergoes 
a compression or expansion; hence if the total duration of the shifted 
waveform is compressed or expanded, approaching the duration of one code 
symbol, energy will be lost in the correlation process. The criterion 
for the use of the narrow-band analysis is given in Appendix B and is 


duplicated here. Hence if 


Ital i 
i Lg NN ; 6-1 
Or alternatively, 
1 Cc * ° 6=-? 
Wi ees (See equation B-14) é 


Z 


Then the narrow-band analysis applies. 





*This result is indicated by A. W. Rihaczekl4, and E. J. Kelly and 
R. P. Wishner.14 
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Equation 6-2 is a statement of the velocity-time-bandwidth-product 
uncertainty relationship. 
B. SYSTEM BLOCK DIAGRAM 

The system block diagram for the wideband case is essentially the 
same as the narrow-band case with the exception of two critical areas. 
A guard interval is necessary between code A and code B to allow for 
time expansion or compression of the two codes due to the doppler effect 
on the envelope of the received signal. Also the delayed replica of the 
transmitted signal must be clocked at the rate corresponding to the 
expected clock rate of the received signal and the carrier frequency must 
be shifted by the doppler frequency. Suppose that the replica signal 
were shifted by the doppler frequency, but not clocked at the expected 
rate; then degradation of the output of the matched filter would be 
expected. In fact, the above statement implies that we are requiring a 
narrow-band system to detect a wideband signal. From the point of view 
of equipment complexity, it 1s more economical just to shift the carrier 
frequency of the transmitted signal to obtain the replica signal. Appendix 
C contains the analysis of the replica signal just being shifted in 
carrier frequency for the wideband case. It is expected that no serious 
degradation of the output of the matched filter, and hence the "envelope 


detector," 


will be evidenced, so long as the time-bandwidth-product-velocity 
uncertainty relationship is not violated. (Equation 6-2) 
C. WIDEBAND AMBIGUITY FUNCTION FOR SINGLE CODES 
In the wideband case the replica signal is obtained by clocking 
the transmitted signal at the expected rate and translating this signal 


by the expected doppler frequency. For convenience the wideband matched 


filter block diagram is given below, 


ie 


yee judi Wd") 


ee ete _ 


Wideband matched filter 





=! (t = tq) 


Figure 6-1 


The replica signal is given by 


$:(D,4- T) = I(D;4-G,) cos w),(D;4- %) ; 6245 


where D. = expected doppler factor 
f;- 
= 4+ 92 
fo 6=3b 
Wai = 


expected doppler radian frequency, 


The received signal is given by 


$(4-44) = D' T[D(A-ty)] cos [C+ (t-ty)] 6-4a 
where D' = target's doppler factor 
f 4? 
-1+— 6-4b 
fo 


Wg! 


target's doppler radian frequency, 
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Then the output from the matched filter is given by 


60 


UY: Ter Wa QD) = ( sit ~t5) $; (Dp b- %) At 


=< 


= o'( T[p'4-4a)] cos [( 0, + wy (4-+4)} © 
TLD, (- FEY] cos [Cuys dys ~ EY ab 6-5a 
Now let us shift the time reference by letting 
ti =t- C4 é 


Thus equation 6-5a becomes 


0 


Mi . O'} T(d4PTLO(4,+t4- ~) cos(wi4 Wy 4, ° 
COS [Cu 4 wy, VO, + 4, - ok] cis Gait 


Then define 


eB 


Che a 
= 6-6 
Ds a ‘ 
Since ty is the variable of integration, it may be replaced by t with 
no loss in generality. Hence the output of the matched filter may be 


written as 


ys = DY TID) TLOL(t-y)] cos (tus) ¢ 
cos (Wri MA-T | db, 6~7a 
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we 


( TWIT (0; 4- TH cos [Coy dysrt 4 (Qs dyed 


¢5 


+2 \ TOOT Lop) cos (aus ay - (as dg gh de | 


-Oo 


6-7b 


The second term in equation 6-/7b will be rejected, provided that 


(Wo' - Wai) & (2W, + Wg"). As in the narrow-band case, the output 


of the matched filter is an ordinate in the T domain. Then, the collator 


collects these ordinates in the € domain. The output of the collator 


is given by 


oS 
Y (re, oda wa’) = D \ Ts) aun (4-7,)\ Cos | (udg- dys) 44 (4 ys) % dt G20 
ov 
Let 
Woy = Wot Wye 6-9a 
oad fe Wy Wt . 6=-9b 


If the Euler identity for the cosine function is used, equation 6-8 becomes 


yltruyz we’) = Of fel" oe Bu To Th Aa-eyje! ue att 6-10 


The output of the “envelope detector" will be given by 


Xqaurysrws’) = Brur Ree fe 5 ( Top Tin ale nat 


Ue 


To the extent that we can ascribe an envelope to the output of the 


summer, then 


Env < We {el 7 ay \ LoIfo, (4-7 yjel eat > = 


60 YH t . 
| 2( Tos) T[0;({-7,] a | 6-11 


—Oo 


In Chapter II, it was proved that for the narrow-band system the output 
of the "envelope detector" is identical to the ambiguity function. For 


the wideband case the ambiguity function will be defined by 


7 00 ° 4 
PG» Way 5 wy’) = | S| T(o# T| 0, (4-71 eye 6-12 


Et 
It is only to the extent that the approximation made in equation 6-11 
holds, that it is possible to give a physical interpretation to the 
wideband ambiguity function. 

Let us now reference the received signal to the replica signal by 


letting 
t; = Dst, 
Thus equation 6-12 becomes 


BY / Dt (a ets 
Ty Mary) = Io, > \s(24,)2U,- Dit) € alt 6-13 
oO 


3 


-_ 


*If the integral in equation 6-11 is a "relatively" slowly varying 
function of ‘%,, as compared with ej%aT. , then this approximation 
is valid. 
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Then define 


D, = 0 = te+ ta’ 
4 D; Sas 4s ) 6-14 
Tice DT, =B =D, 6-15 
_ Wd _ /tdy’- a 
Us, = > 7 [a= dt ) ie 6-16 
£ Wo + Was 


= (D-L)® , 6~16b 


Therefore the wideband ambiguity function becomes 


Y Da ( jt 
(7, W,.) = 2s \ T(Dyti TA ODe? th 6-17 


-OoO e 





Note that either equation 6-17 or equation 6-13 may be chosen as the 
definition of the wideband ambiguity function. 

One of the properties of the narrow-band ambiguity function is that 
it is symmetrical in range deviation, T , and doppler, W,. However, 
the wideband single code ambiguity function seems to exhibit no simple 
symmetry property. Thus, we will obtain equations for the wideband 
single code ambiguity function, explicit for a limited case,to indicate 
the method that is involved. 

The wideband single code ambiguity function will be calculated for 
w, 7 0, T>O. These equations will be limited to the case of the 
received signal's envelope being compressed or expanded by not more than 


the duration of one code symbol of the replica signal. Appendix E 


id 


contains the analysis of the Limits of integration and the code values 

in the interval of integration for the wideband, single code ampiguity 

function. We have divided the ambiguity function into two intervals of 
range deviation. 


1. For range deviations, 'T , in the interval 


, / To 
(m-3) 2 = = ey ~ Nad (mye ; MM eA,Z,o0N, 





1 


we have 
Gap’ 
4 | N- 3 
! « DRC Ss {at 
Cz, w,) = > ) X;Xiap | edit 
ee TAi-W) To, 
Ne(pas) Ta Noi cae 
ats : Ky spas e! dt 6-18 
- Gepyr’ > 
where N = length of the code >) 
xX. = Bal depending on the ith value of the code symbol 
Dees veplice signal's doppler factor = 1 + a 
O 
t 
D' = received signal's doppler factor = 1 + fd 
O 
t f,' 
DA i dD. = aod 
D o + fai # 


Wa = (Dy, — 1) Wo 


T = transmitted signal's duration of a code symbol 


i 
Tp; = replica signal's duration of a code symbol ata 

at 
Ty’ = received signal's duration of a code symbol = <> 


Ty) 
Dy 
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AT 


magnitude of the difference in the duration of a code symbol 











a _ Tt 
i Tye rT 
p = smallest integer > Ire | 
6-19 
j = ‘|= (denotes a complex quantity) (cont.) 
Z : ; . |2 
Q = normalization constant ES 
Then if the integration is performed, equation 6-18 becomes 
t 2 2 
rane W),) = DLO — KKy, Gaeyrt G-a toi - TI » 
mi 2 
wal lisp + G-0Tos 
pin WelGapT-G-aTp-T) {SA ek 
Zi | S 
W1G pT - G-Tp.-74 
Zz 
N- (pat) . 
+2, x Lica ity p7jsin ale wd 
= KiKi spas LCr Mpj- Gapyt to Ze 
a Wal Tes tp:- aT d 
Ar 
© 1s Li Tp; + Ga pT j 
e! 6-20 
2. For range deviations, @ , in the interval 
: f 
oa yt a ee : ae 
cs N+i4(m Ve Ss ee m= 1,2,°°,N; 
we have 
' j 
4k) 
2 | n-k-k mond: 
! OAD 8 
Bawy = MAN TD xx) |e “at 
- T+ LY Toy 
n-§- (eed) T+ iT; nl 
— ; A 
ed Jerid c dé 
a Gak)yT’ 


ie, 


! 
Cn-8aLyT 
-Wyt 


a Ang Xn gar \ edt 


(n-Ryr ! 
ne 4 G-kK)TG; aon 
{ \ 4 
ata, Xe (\ e dé 
aa (ies) 
N2 Ga2yT” es 
| é 
+2 Miter Kive \ edt = 
4= f- 
: T+Ci-k) Tp, ’ 


where the definitions in equation 6-19 apply and 
Bs 


AT 


n 


_ Ipi 
ny = Me 


q = smallest integer > sr 


k= 
tt 


q modulo (ny) 


nn 
i 


smallest integer > T/Tou| 





If the integrals in equation 6-21 are evaluated, then the ambiguity 


function becomes 


n-0-k 
Pia wy = 0px? | 2, XX, skort Gy y-7] ° 
a 


° 
4= 





ain, WalGiskor! G-sT-7] : WalGskyT 4 G-)T0a 


WL Gk T- (-1T; -T) 
vA 


n-0-(K41) 
a um PGP rt | T+ ilo: _ Gakyr' | S 
4=1 2, 
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— Ore —-GaldT’] ey LTT Gao] 
sin | me @ J 2) | ee 
Wo [e+ ATp-G4kIT J 
Zz 


, / 
WaT > Wall@n-2044 )T-c] 
‘Sin AL 
a Ket X ae a wea a} 


xT 
Zz 
N-Z 
+ ee eee Catt, SU90m ls 
1= n-Q 7. 


a (t+ Gk -<GaT | wl G-k) To: Gar J 
Sin ae a | e! 
WalT+lUi-k) Ds- G4aT | 
Z. 





4= f- 


N-2 
° wee Sig 
+ Z, Mi est oe (442 \ =a Eire) o 


Ww) Geer - (iak) "pb; = T | © Wa [(Cia2)To (4-)Tp; 
ey e 6-22 


Ws (DT ~ Gab T,-T 
pF 


WIDEBAND AMBIGUITY FUNCTION FOR COMPLEMENTARY CODES 


As stated in Chapter IV-C, to obtain the ambiguity function for a 


complementary code, given the ambiguity function for single codes, one 


merely replaces the signed product for the single codes by the signed 


products for the complementary codes. 


Consider range deviation, ' , in the interval 


a L ez nee Tos 
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Then if the signed products for the single codes in equation 6-20 are 
replaced by the signed products for the complementary codes, the ambi- 


guity function for the complementary codes will be given by 
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where ix = =, 


2. Now consider range deviation, T, in the interval 
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then if the signed products in equation 6-22 are replaced by the signed 
products for the complementary codes, the ambiguity function for the 


complementary codes will be given by 
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CONCLUSIONS 


A. SUMMARY 

The research presented in this thesis had two primary objectives. 
The first objective was to formulate the ambiguity function of the comple- 
mentary codes for the narrow-band case. The second objective was to 
formulate the ambiguity function of the complementary codes for the 
wideband case, We shall summarize the results of these investigations, 
but not necessarily in the order presented in the thesis. 

1. NARROW=-BAND ANALYSIS 

In the narrow—band case, Siebert's ambiguity function was extended 
to include complementary codes. The ambiguity function of complementary 
codes was defined so that its autocorrelation function contained no range 
sidelobes. From the definition we determined that the two codes should 
be combined at the outputs of the matched filters to effect the necessary 
cancellation of range sidelobes. During the course of research into the 
ambiguity function of complementary codes, we discovered that, once the 
ambiguity function for a single codes has been derived, then to derive 
the ambiguity function for a complementary code is a simple task -- replace 
the signed products of the single code with the sum of the signed products 


of sequences A and B, 1-e., 


replace (X;X5) by (XA; XA; +X 4 XB; ) . 


In the derivation of the ambiguity function for single codes, it 
was proved that the ambiguity function was identical to the output of 
the "envelope detector." (The “envelope detector" refers to the device 


in the system that takes the envelope of the signal, but in the range 
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or TT domain. This use of the term differs from conventional terminology 
in that, usually, envelope detection is done in the time domain.) Although 
this result is not unique, it helps to clarify some of the misinterpre- 
tations of the ambiguity function. It is a commen misconcepticr to state 
that the output of the matched filter is identical to the ambipuity 
function, 
One result obtained in the analysis of the ambiguity function of 
complementary codes is that in the doppler domain, the zero crossings 
of the ambiguity function move closer to the origin as the length of the 
code increases, if the duration of a code symbol is held constant. The 
implication of this result 1s that longer codes permit better velocity 
dierrimination. However, this process of increasing the code length to 
obtain better velocity discrimination cannot be continued indefinitely 
without incurring a penalty. The penalty is that if the velocity-time- 
_bandwidth-product uncertainty relationship is violated, then distortions 
of the envelope of the received signal must be taken into account. This 
amounts to having to clock the replica Signal at the expected rate and 
hence increasing the complexity of the system. 
2. WIDEBAND ANALYSIS 
In the wideband analysis, it can be shown that the doppler effect 
is three-fold, 
a. The carrier frequency is shifted by the deppler frequency. 
b. The envelcpe of the transmitted signal is compressed or 
expanded according to whether the target is closing or 
opening with respect to the receiver, 
c. The received energy differs from the transmitted energy 
by the same factor as the received and transmitted carrier 


frequencies. 
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Although all three effects can be proved from the special theory of 
relativity, only the last two were proved in this manner, The distortion 
of the envelope of the received signal was also proved as a logical con- 
sequence of the application of Fourier analysis, and the doppler effect 
on the carrier frequency. A. W. Rihaczekl2 proved, by considering a 
photon, that the received and transmitted energies differ. Our approach 
to the problem was to make a Lorentz transformation of the field tensor. 


It was proved that for binary phase-reversal modulated signals, if 
Cc 
NTW|vy] <€ > 


then the narrow-band analysis applies. 

The wideband ambiguity function was derived here; however this 
result had been previously accomplished by A. W. Rihazekl4, and Kelly 
and Wishner!4, We merely clarified certain aspects of the analysis. 
Then the equations for the single code and complementary code wideband 
ambiguity functions were derived. These equations are limited to the 
case of signal's envelope being altered by not more than the duration 
of one code symbol. They are, indeed, formidable, and a computer is 
required to obtain the ambiguity diagrams. Due to time limitations, the 
wideband ambiguity diagrams were not completed. 

B. AREAS OF POTENTIAL STUDY 

Many areas of future research follow directly from this thesis. 

A logical extension of this thesis would be to construct the wideband 
ambiguity diagrams for complementary codes and single codes (Barker codes, 
truncated pseudo-noise codes, etc.) Also the region of validity for the 
narrow-band analysis should be more precisely defined. We found that, 

in the preliminary analysis of the computer runs for the wideband ambi- 


guity function of complementary codes, the wideband ambiguity function 
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for complementarv codes began to deviate substantially from the narrow~ 


band ambiguity function at approximately 


ntw |v, = 0.1 ¢/2 
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APPENDIX A 
ENERGY OF THE RECEIVED SIGNAL 
DERIVED FROM CONSIDERATIONS 
OF THE SPECIAL THEORY OF RELATIVITY 
Consider the reference frame, S, with the transmitter located at 
the origin of the coordinate system. Assume that the target is a perfect 
reflector, which is moving at a velocity, v,, relative to the transmitter, 
The target is situated at the origin of the reference system, S'. A 
Suitable set of axes is chosen so that the x direction of both of the 
systems, S and S', is collinear with the relative velocity vector, v_. 
figure A-l reflects a suitable choice of axes. Note that in Figure A-l, 
the convention for positive velocity is away from the transmitter and is 
a standard convention for most physics textbooks. Although this con- 
vention is contrary to the one chosen in the main text of this thesis, 
no confusion should arise here, 
The transmitter emits a plane wave of finite duration, NT. Let field 
quantities as seen in the transmitter reference frame be Fy and H,. 
Next, we require the field quantities as seen by the target in 


reference system, S'. These quantities can be obtained by application 


of the Lorentz transformations of the field=strength tensorl-, They are 


given by 
1. te a = ' — = 
La = 0 Hy H, = 0 
7 Ff = ¢ _ ' = ° = om 
Ey BCE, - BH) Ke Oe) = 0 A-1 
Ral ee : aes t = 08 7 
E, ¥(E, + BHO i Bz - BE) , 

where 
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Reference frames for the target and the transmitter 


Figure A-1 





The plane wave after reflection 


Figure A-2 
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The plane wave after reflection 
as viewed by the transmitter 


Figure A-3 


90 


Then we want the reflected wave's field components as viewed by 
the target. These field components can be easily obtained by applying 
the boundary conditions for a perfect reflector. Therefore, the field 


quantities as seen by the target are given by 


ao He a= 0 

ee o ee hy! = 0 ead 
ee) H'' = 4H ' , 

Zz Z Zz 


Finally, we require the field quantities as viewed by the transmitter. 
The reflected wave's components can be easily obtained by applying the 
transformation rules in equation 1, and noting that the velocity component 
is now in the negative x direction; hence we must replace vy by -v_. 
Let us denote the reflected field quantities by a triple prime notation. 


Thus 


2a = 0 oe = Q 

es = ane ae Qu,'") ee = () A#W5 
a ter — tt ~ tt 

Bt = 0 He = BC + BELT), 


If equation A-4 is substituted into equation A-5, 
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1 = a oe A-6 
Then substituting equation A-1 into equation A-6 
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The solutions for Maxwell's equations are such that 
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since ns ee = 1, if the field quantities are in Gaussian units. 


Therefore, equation A-/ becomes 


Ey''' = [-¥*a + 87) + 2x28] E, 


yi 
= -¥2(1 - @)°E 
y » A#W9 
aan 2 2 2 
Hep Behe Oa. —) 2am 
= Sa 
= O°(1 - ®)7H, . A-10 
The power density can be obtained by Poynting's Theorem, and is 
given by 


p = 5 Rei Ex H*{ 5 A-11 


Thus the power density of the reflected wave as seen by the transmitter is 
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where Ey is the amplitude of Ey 
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The power density of the transmitted wave is 


Pl = lel i 


The energy density of the wave can be determined by integrating 
the power density over the duration of the signal. Note that the power 
density is uniform over the planar surface. Thus the energy density of 


the transmitted wave is 


GO 
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And the energy density of the reflected wave is given by 


OO 
val 
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where N = length of the code 
T''' = duration of a code symbol for the received signal 


T = duration of a code symbol for the transmitted signal, 


During the analysis presented in Chapter V and Appendix D, it is shown 


that the received signal's duration of a code symbol is altered by the 


doppler effect so that 


peers == 


DpH 


where 


D = vz Vise 
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Thus equation A-15 becomes 
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Then the ratio of the received signal's energy density to the transmitted 


signal's energy density is given by 
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i DEP Te ND Te aaa 
CRITERION FOR THE USE OF THE 
NARROW=BAND APPROXIMATION 
(VELOCITY-TINE-BANDWIDTH-PRODUCT UNCERTAINTY RELATIONSHIP ) 

In this appendix we shall derive, from rather elementary consider- 
ations, the region of validity for the narrow-band approximation. This 
approximation is made in Chapters II and III and assumes that doppler 
has negligible effect on the envelope of the signal. Chapter V contains 
the derivation of the time waveform of the doppler shifted signal for 
the wideband case. (For purposes of clarification, the terms, doppler 
shifted signal and the received signal, are used synonymously.) We con- 
cluded that the effect of doppler was threefold: 

1,. The duration, Tp, of a code symbol for the doppler shifted 

waveform was altered so that 


Tp = TAs B=-1 


2. The carrier frequency was shifted by the doppler frequency, 
3. The received signal's energy density differed from the trans- 
mitted signal's energy density by the doppler factor. 

If the doppler shifted signal's duration of a code symbol is altered by 
doppler so that the total duration of the signal 1s expanded or compressed 
by the duration of a code symbol of the transmitted signal, we expect 
degradation of the output of the matched filter, In particular, if the 
code is completely random, then for zero range deviation the output of 
the matched filter will be zero. Hence we shall use as our criterion of 
the narrow-band approximation, that the total duration of the received 
sipnal must be compressed or expanded by much less than the duration of 


a code symbol of the transmitted signal, T. Let us examine two cases, 
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1. Case I: The total duration of the received is compressed by 
much less than the duration of a code symbol of the transmitted signal. 
(This statement implies that positive doppler is present.) 

Figure B-l indicates the envelopes of the transmitted signal, I(t), 
and the received signal, I(Dt), for some arbitrary code. kence for the 
condition that the total duration of the received signal is compressed 
by much less than the duration of a code symbol of the transmitted signal, 


the following inequality applies: 

NT, YD (N-1)T. P=-2 
If the above inequality is rearranged, 

TE econ (N = 1) E-3 
From equation Bel 
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Note that the time-bandwidth product is equal to }:.* 


Hence 


NTW v_ We C/Z2 B-6 


2. Case II: The total duration of the received signal is expanded 
by much less than the duration of a code symbol of the transmitted signal. 
(This statement implies that negative doppler is present.) 

Figure B-2 indicates the envelopes of the transmitted signal and 
the received signal for some arbitrary code. Hence for the condition 
that the total duration of the received signal is expanded by much less 
than the duration of a code symbol of the transmitted signal, the following 


inequality applies: 
NT D972 (N- 1)Tp. B-7 
The above inequality may be rearranged so that 


T/Tp >» (-1/N. B=8 


Note that the doppler frequency is negative in this case. Hence the 
doppler factor may be written as 


\£q\ 


£ 


bz=1-.- 





O 


*We have defined bandwidth as the video bandwidth to the first null 
of the rectangular pulse spectrum, i.e., W = 1/T. This is approximately 
equal to the "effective'’ RF bandwidth. Thus NIW = N. The time-band- 
width product, as used by Torres, was defined in such a fashion (using 
standard deviation in both time and frequency domains) as to be a complete 
and precise description for a Gaussian pulse. If rectangular or other 
shaped pulses are used, a numerical factor of the order of magnitude of 
unity may be introduced into his precise inequality. Hence there exists 
an apparently discrepancy of a factor of 2 between Torres’ inequality 
and our inequality B-6. 
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If equation B-9 is substituted into the inequality P-&, then 
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Alternatively, the doppler factor may be given by 


l \vri\ 
D = 2 oe 


ie, B-11 
c 
Hence, substituting equation B-1l into the inequality B-8, then 


ae Z i 
—— S3F 7-7. B-12 


If the length of the code is much greater than unity 


\vyz\ yds 
C « 2N e B-13 
Then, we note that the time=-bandwidth product is equal to the length 


of the code, N. Hence 
c 
NTW wy x 7° B-14 


Thus if the time-bandwidth product multiplied by the radial component 

of relative velocity is much less than one-half the velocity of the wave 

in the medium, then the narrow-band assumption applies. Inequality B-14 

is a statement of the velocity-time-bandwidth-product uncertainty relation- 


ship. 
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APPEND EC 
ERFECT OF DOPPLER SMISMALCH SsDUr Om [hija lich 
SIGNAL ONLY BEING TRANSLATED IN FREQUENCY, ON 
THE OUTPUT OF THE "ENVELOPE DETECTOR" 
1. BACKGROUND 

For the wideband case, the system block diagram contained range 
banks, such that in one range bank, the received signal is correlated 
with an exact replica of itself. Note that the replica signal has the 
same carrier frequency and clock period as the received signal. llence, 
in generating this signal from the transmitted signal, the system must 
contain a clocking device, such as a shift register, as well as a fre- 
quency translator, Naturally the clocking device adds to the complexity 
and hence the cost of the system. Suppose that in trying to reduce the 
complexity of this system, we merely translate the frequency of the 
transmitted signal by the doppler frequency. How is the output of the 
"envelope detector" affected? If the replica signal is obtained by a 
mere frequency translation, we are, in effect, requiring a narrow-band 
system to detect a wideband signal, Thus we expect no serious degradation 
of the output of the envelope detector, so long as the time-bandwidth- 
product uncertainty relationship (equation 6-2) is not violated. 

2. CUTPUT OF THE ENVELOPE DETECTOR (SINGLE CODE SYSTEM) 

We will investigate the single code system prior to the complementary 
coded system, since once the equations for the single code system are 
obtained, it is a simple task to obtain the equations for the complementary 
coded system. Figures 2-2a and 2-2b will be used as the system block 
diagram, since our receiver is a narrow-band system, The received signal, 
s'(t), will be given by equations 5-10a and 5-10b. The replica signal, 


S(t), is simply the transmitted signal translated by the doppler frequency} 


LOL 


hence it is identical to the narrow-band received signal. Then s(t) 


is given by equation 2-4, Thus 


s(t) = D'T(O'4) cos(wW+wat , C-ae 


S y(t) = T(4) cos (We+ DIE C-1b 
Then the output of the mismatched filter in Figure 2-2b is given by 


/ 60 
Yay (> Wy: oY) : ” cos wt § (oe) 14-7) dt 
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‘aa TOA) TU-) cos] AC, +0916 - (4 dye] db, C-2 


The second term in equation C-2 will be rejected by the integrator. 


Then the mismatched output of the "envelope detector" will be given by 


XK yb™ mol & | (roe) TU-2) dt, | : C=3 
UL = wa’ 


From equation 6-17 we see that the mismatched output from the "envelope 


detector" is related to the wide-band ambiguity function by 


Ky Se = (%0,) C-4 


Vrel 
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Recall from equation 6-14 that 
Dd, = eae 


Then the differential doppler factor, in this case, is 


Dae C=5 


Next, let us obtain the mismatched output of the "envelope detector" 
for the case of the replica signal's doppler frequency matched to the 
received signal's doppler frequency. ‘This can be obtained by setting 
ve ft, and Ww, = QO, in equations 6-20 and 6-22. 


1. For range deviations, T , in the interval 
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where N = length of the code 


eee 1 depending on the ith value of the code symbol 
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D' = received signal's doppler factor 
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T = transmitted signal's duration of a code symbol 


Tp' = received signal's duration of a code symbol 
= T/D' 
AT = magnitude of the difference in the duration of a code C-7 
(come. ) 
symbol = |T - T,'| 
p = smallest integer >\ T/-\ 2 
2. For range deviations, tT, in the interval 
a ! 
oe Nee 2 Glo me 142,068 yk 
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where the definitions in equation C=-7 apply and 
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q = smallest integer > \ "c/a | 
C=9 


smallest integer 2 \"/+\ ; 


3. OUTPUT OF THE ENVELOPE DETECTOR (COMPLEMENTARY CODED SYSTEM) 

As indicated in Chapter 4-C, it is easy to obtain the equations for 
the complementary coded system once the equations for the single code 
System are obtained by simply replacing the signed products of the single 


codes by the signed product for the complementary codes. Hence if 


XiXy is replaced by XAGXAZ + XB ZXBy 


in equations C-6 and C~8; then the resulting equations will be the mis-~ 
matched output of the envelope detector for the complementary coded 


System. 


LOS 


APPENDIX D 
DURATION OF A CODE SYMBOL DERIVED FROM CONSIDERATIONS 
OF THE SPECIAL THEORY OF RELATIVITY 

In Chapter V it was proved that the transmitted signal's duration 
of a code symbol, T, differed from the received signal's duration of a 
code symbol, T', by the inverse of the doppler factor, D. This result 
was obtained from considerations of the doppler effect on the carrier 
frequency and also by the use of the Fourier transform. 

In this Appendix we will prove the same result, but from consider- 
ations of the special theory of relativity. 

Consider the reference frame, S, with the transmitter located at 
the origin of the coordinate system. Next, consider a moving frame, S', 
with the target located at the origin. A suitable set of axes is chosen 
so that the x direction of both S and S' is collinear with the velocity 


vector, v At some initial time, 1.e., t = 0, we shall assume that the 


re 
two frames are coincident. At t = 0, the transmitter emits a plane wave 
pulse of duration, T. The pulse travels at the speed of light. Then we 
can view the pulse as extending in space so that its trailing edge has 

an x coordinate of -cT. For t > 0, the moving frame travels in the +x 


direction with a velocity v Note that the convention for positive 


ee 

velocity is away from the transmitter and is a standard convention for 

most physics textbooks. Although this convention is contrary to the one 

chosen in the main text of this thesis, no confusion should arise here. 
At time t = T, the trailing edge of the pulse reaches the origin 

of S, and the leading edge of the pulse is now at x = cT. Now let us 


consider the spacial distribution of the pulse in the moving frame. In 


the moving frame the x coordinate, x', and the time in this frame, t', 
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can be derived from their counterparts in the reference frame from the 


Lorentz transformations. ~° Thus 


x' = Y¥(x - vit) D-1 
t' = ¥(t — sew), D=2 
c2 


where 


vo 2 -1/2 
ae ale ) : 
Let us denote, as in Figure D-l, the trailing edge of the pulse as event l. 


For event 1, the x coordinate of the trailing edge of the pulse in the 


moving frame is given by 


4 
Ui 


: S(x - v.t) 


rt 
il 
r4 


“ oven ‘ D=3 


For event 1, the time elapsed in the moving frame 1s given by 


Vv 
ta $(t - 7 *) x 


il 
o 


ST D-4 


Since the pulse is traveling at the speed of light, the total time for 
the trailing edge of the pulse to pass the observer at the origin of the 


moving frame is given by 


Oy: 





Reference frames for the 
transmitter and the target 


Figure D-~1 
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Reference frames for the 
target and the receiver 


Figure D-2 
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Equation D-5 is a statement of the expansion or compression of the 
duration of a pulse for the case of one-way doppler. 

Now we would like to consider the duration of the pulse which is 
reflected from the target, as seen by an observer (the receiver) at the 
origin of the reference frame, S. To avoid confusion let us denote this 
reference frame as S''s; however, bear in mind that we are assuming that 
the transmitter and receiver are at the same physical location. We can 
view the reflected wave as being retransmitted from the target. 

Consider S', the target's frame, and S'', the receiver's frame, to 
be coincident at t’' = 0. From S', the receiver seems to be moving at a 
velocity v, in the negative x direction. Then at t" = 0, the leading 
edge of the pulse is at the origin of both frames. At t' = T', the 
trailing edge of the pulse is at the origin of the frame, S', Let us 
denote this as event 2. To obtain the x coordinate of the trailing edge 
of the pulse as seen by the observer on S'', we utilize the Lorentz trans- 
formations as given by equations D=-1 and D=-2, but remembering to replace 


vy, by -v,- Thus for event Z 
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For event 2. the time elapsed in S‘' is 


t, = 6(t'- Ex’) |, 


X=0 
4=7 


Sai D=-7 


Then the total time required for the trailing edge of the pulse to reach 


the observer on S'' is given by 
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Then if equation D-5 is substituted into equation D-8 
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= T/D , D-9 
where 
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The factor which compresses or expands the duration of the pulse is 
identical to the one obtained by the Fourier transform analysis done in 
Chapter V, with the exception of the difference in sign of the velocity 
due to the change in sign convention adapted in this Appendix. In our 
initial assumptions, we assumed that the duration of the transritted 
pulse was T. It 1s easy to see that if a series of N pulses, each of 
duration T, were transmitted, then the compression or expansion factor 
would be identical to the one in equation D-9. Hence we have proved that 


the doppler effect alters the duration of a code symbol by the factor D. 
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APPENDIX E 
DETERMINATION OF THE LIMITS OF INTEGRATION 
FOR THE WIDEBAND SINGLE CODE AMBIGUITY FUNCTION 

In this Appendix we would like to determine the limits of integration 
and their corresponding code values in the interval of integration for 
the wideband single code ambiguity function. In the formulation of the 
equations for the wideband single code ambiguity function, there are 
four quadrants to consider. We will consider only the first quadrant, 
ioe, Ww, 70, 70, and will limit the validity of the equations to 
the case of the received signal's total pulse duration being expanded 
or compressed by not more than the duration of one code symbol of the 
replica signal. 

By equation 6-17, the wideband single code ambiguity function is 


defined as 


euch 
Y'(7, d,) = \%\to,oTa-2) ade . E-1 
-00 


Hence to determine the limits of integration and the code values in the 
interval of integration, we must analyze I(D, t) and I(t =). “Doe 
is the envelope of the received signal (time referenced to the replica 


signal), and is given by 
N ( 
L(t) = Z| ALCv4-i79)- Mo.t-GaayTo i, E-2 
el 


Then define 


Th = Tp /, r~3 


jee? 


Therefore, equation E-2 becomes 
N ; 
TD,k) = ADK, 4 Alb-it’) - Me- Gay i 5 ie 
c= 1 
The envelope of the replica signal is given by 


N 
mic) =I 1 fllt-- itp; ) - {\{t-e- G+ Tot ; bes 
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Figure E-l is a set of waveforms of the received and replica signals' 
envelopes. We have chosen, in Figure E-l, a code of length 6, with the 
code values of eee =I, atl, +1. -i% +1} as an example. From this example 
we will obtain an equation that will apply, in general, to codes of any 
length. Note that the amplitude of the received signal's envelope in 
Figure E-l is less than the amplitude of the replica signal's envelope. 
This is done deliberately so that both waveforms may be seen in detail; 
however, bear in mind that both envelopes are, in actuality, of the same 
amplitude. 

For positive doppler we recall that the received signal's waveform 
will be compressed. Figure E-l is a series of waveforms of I(@D, t) and 
H(t), invwhich we hold [(D, t) fixed and shitey((t.— {) to tieeaeene, 
corresponding to positive range deviation, T . In doing so, we see that 
the limits of integration (as a general rule) change for 'T in increments 


of AT, where 


AT = ae = 7) ° E-5 
Next, we would like to express the interval of integration (and hence 


the limits of integration) in the form of a series. Let us now look for 
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symmetries in Figure E-1. Note that for 04'T4 4 AT, the limits of 
integration over each interval, in addition to the code values of the 
received and replica signals, remain constant. We also see that this 
occurs again for 9 AT4 T <13 AT. However, we have lost the duration 
of one code symbol of the replica signal. Similarly, for 19ATZ@T< 23AT, 
we see that we have lost the duration of two code symbols. From this 
analysis, we can obtain a general formula to express the limits of inte- 
gration and the code values in the interval of integration. Hence for T 


in the agnterval 





/ 
(m5 £ ef ‘ -N+1- (m+a) “Be 5 m= 1,2,°°°K, 
where N = length of the code 
= +1 depending on the ith value of the code symbol 
D, = replica signal's doppler factor he 
D' = received signal's doppler factor 
D, = differential doppler factor = D'/Dj 
Boas Da 
7 = Tgp 
AT=\T, -T'| , 
The wideband single code ambiguity function becomes 
GapyT 
%'Ce) Ws) = O00 | Sx, See eas 
= w+ G-2) Tp, 
N-(p+1) w+ iToe is 
+> X; LA fapar \ e * dt ,) 7 
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in which the definitions in equation F=-6 apply and 


p = smallest integer > \ fa! 
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normalization constant = \ 2 
wd - 0 | 
Wr = ( Aan~ ye = (D- 108, . so 
Udo s VAL 
To complete the analysis, we would like to obtain the code values 
and limits of integration for the intervals, 4 AT4T 4 BAT, 14 ATLTCIBAT, 


and so on. As a general rule, the interval of range deviation becomes 
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We see that we may express the ambiguity function for this interval as 


the sum of five terms, 1.@o, 
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where the definitions in equations E-6 and E-8 apply and 


2 
AT 
‘TE ° 
n= 
AT 
q = smallest integer E-10 
1 = q modulo (ny, ) 
k = smallest integer > \ Gia | 
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PROPERTIES OF COMPLEMENTARY SEQUENCES 3, 15 


This Appendix is a compilation of some of the properties of comple- 
mentary sequences, At the end of this appendix is a list of the known 
kernels, of which there are four. 

1, Definition: A pair of binary sequences of equal length, with 
the number of like pairs of one sequence, A, equal to the number of 
unlike pairs of the other sequence, B, for each possible spacing except 
at zero spacing, where the pairs of both sequences are like, is said to 
be a complementary sequence. 

2. Given sequences A and B each of length n. Then complementing 
the code, 1.e., complementing each element in the code, of A or B or 
both, results in a complementary code. 

3. Definition: Time inversing of a code is to interchange the lst 
and last bit, the second and next to the last bit, and so on. 

4, Time inversing the A code, B code or both results in a comple- 
mentary sequence. 

5. Definition: The altering of a code is a transformation where 
every other bit of both codes is complemented. 

6. The result of altering a complementary sequence is again a 
complementary sequence. 

7. Definition: A kernel is a basic length code which cannot be 
decomposed into shorter length codes by an inversion of the standard 
generating methods, 

8, Definition: Complementary sequences which are not kernels are 


called composite complementary sequences. 


il? 


9. If A = ajscccevoesay 
B= Dy scccccces by 
is a complementary Sequence pair, then 


é = Aprecccccesd sD] yocveccoyhy 


— x 
D = Apssescccesand. poocccoeg Dy, 
form a complementary Sequence pair. 


Loe if A = Ajseeccccoesd, 


B= Dypececccosb, 


1s a complementary sequence pair, then 
C = Ayedy sagsbdorceccsooeeosansbdy 
D = ais Biiapulg 1s eee eee onRen 
form a complementary Sequence pair. 
11. The following is a list of the known kernels: 


ae length 2 


ieee 
eee 


A 
B 


b. length 10 


A, = {-1,41,41,-1,41,-1,41,+1,+41,-1} 
By = 1-1,+1,+1,+1,+1,+1,+1,-1,-1,+1¢ 


Aj = (41,-1,+1,-1,+1,+1,+1,+1,-1,-1| 
By = {+1,+1,+1,+1,-1,+1,+1,-1,-1,+]} 


ce length 26 


A=. {+1,=1,£1,+1,=1,-1,42,-1,-1 9-141 
+1,-1,-1,-1,-1,+1,+1,-1,-1,-1,+1,-1,+1 } 
B= {-1,+1,-1,-1,+1,+1,-1,+1,+1,+1,+1,-1,-1, 
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*Bj; indicates the complement of bj 
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